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Abstract 

pi. . Abstract. The major sources causing deterioration of optical quahty in extremely large optical telescopes are 

f— ^ ' misadjustments of the mirrors, deformations of monolithic mirrors, and misalignments of segments in segmented 

^S) i mirrors. For active optics corrections, all three errors, which can partially compensate each other, are measured 

^ ' simultaneously. It is therefore of interest to understand the similarities and differences between the three 

^ . corresponding types of modes which describe these errors. The first two types are best represented by Zernike 

polynomials and elastic modes respectively, both of them being continuous and smooth functions. The segment 

\^ . misaligment modes, which are derived by singular value decomposition, are by their nature not smooth and 

in general discontinuous. However, for mirrors with a large number of segments, the lowest modes become 

effectively both smooth and continuous. This paper derives analytical expressions for these modes, using 

C/3 ' differential operators and their adjoints, for the limit case of infinitesimally small segments. For segmented 

mirrors with approximately 1000 segments, it is shown that these modes agree well with the corresponding 

lowest singular value decomposition modes. Furthermore, the analytical expressions reveal the nature of the 

Q '. segment misalignment modes and allow for a detailed comparison with the elastic modes of monolithic mirrors. 

• ' Some mathematical features emerge as identical in the two cases. 

O 

^ ■ 1 Introduction 

The first large optical telescope with a segmented mirror was the ten-meter Keck Telescope. Its primary 
mirror consisted of 36 hexagonal segments. All currently envisaged extremely large optical telescopes will have 
segmented primary mirrors, but with a much larger number of hexagonal segments, of the order of one thousand 
C^^ I or more. Assuming that the shapes of the segments are perfect at all times, the only source for wavefront errors 

^^ ' generated by a segmented mirror would be misalignments of the segments. These can be measured either 

optically or by edge sensors at intersegment edges. Both types of measurements can, in principle, detect the 
^^ ' relative piston and the relative tilt of adjacent segments. The misalignments are then corrected by appropriate 

(^ I vertical position changes of actuators which can modify piston, tip and tilt of each segment [1]. 

fj ' If the signals given by the measurements were free of noise and the actuators were perfect, the misalignments 

•^ , of the segments could be perfectly corrected. But, in reality, both signal noise and correction noise will give 

^i^' rise to residual misalignments. Both the signals and the actuator positions can be represented as vectors with 

f~i I suitable sets of base vectors. If the errors are introduced by the actuators themselves, one can choose an 

Mh' arbitrary orthogonal set of base vectors. However, if the errors are due to the signals, an optimum modal 

correction requires the possibility to expand the signal into a set of well distinguishable, that is orthogonal 
signal modes, which correspond to well distinguishable, that is orthogonal, actuator modes. Such modes are 
, . , the optimum alignment modes. For mathematical reasons, also described in section 2, they were called normal 

j^ ■ modes in reference [1]. These modes can be calculated by Singular Value Decomposition (SVD). In the rest of 

the paper we assume that the actuators are perfect and noise is only introduced by the signals. 

For segmented mirrors with a thousand segments or more the modes with the lowest so called singular values 
are rather smooth. They are also the ones which appear statistically with the largest amplitudes in an expansion 
of the misalignment error due to random signal noise. In modern large active optical telescopes [2] containing 
both segmented and flexible monolithic mirrors the lowest normal modes could largely be compensated or 
corrected by other elements in the telescope optics. This shows the need to compare the normal modes with 
other sets of modes used in large telescopes, for example Zernike polynomials or elastic modes of monolithic 
mirrors [3]. 
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In the limit of an infinite number of segments over a finite area the lowest normal modes should be continuous 
smooth functions. If the relationship between a position function and the corresponding signal function is 
expressed by a continuous operator, such smooth functions can be calculated from differential equations and 
appropriate boundary conditions. It should therefore be possible to derive analytical expressions for the lowest 
normal modes in the limit of infinitesimally small segments. 

Section 2 outlines the method to obtain the analytical solutions. Before treating the case of two-dimensional 
segmentation the method is applied in section 3 to two examples of one-dimensional line segments. In section 4 
the method is used to calculate analytical expressions of the normal modes of a mirror with hexagonal segments 
and section 5 compares the analytical solutions with the SVD solutions of a highly segmented mirror. Section 6 
discusses the relationship to modes calculated with variational methods to minimise the r.m.s. of the signal 
vector for a given r.m.s. of the actuator vector and section 7 shows the relationship to elastic minimum energy 
modes of a circular plate. 



2 Method to obtain analytical solutions 

A given set of actuator positions is described by a vector v in a riact-dimensional vector space A, and a given 
set of signals by a vector u in a risig-dimensional vector space S. Without noise effects, that is if the signals 
depend only on the actuator positions, the actuator vectors are related to the signal vectors by a linear mapping 
A from A to S, represented as a two-dimensional rigig x riact matrix A : 

Av = u (1) 

In this case a misalignment could be perfectly corrected by actuator changes. In reality, noise in the signals is 
not negligible and only those signal vectors which are in the image space of A under the mapping A can be fully 
corrected. An optimum modal correction of a misalignment of the segments requires the possibility to expand 
the signal in an orthogonal set {u^} of signal modes which are related by the mapping A to an orthogonal set 
{vi} of actuator modes. Such sets are generated by the Singular Value Decomposition of A : 

A = USV" (2) 

U is a matrix of dimension risig x riact and S and V^, where the superscript H denotes the adjoint matrix, are 
square matrices of dimension riact x riact- The matrix S is diagonal and contains riact so called singular values 
Csvdji- The riact column vectors of U and row vectors of V form orthonormal sets {ui} and {v^}, respectively, 
that is, with <, > denoting the inner vector product, 

< Ui, Uj >=< Vi, Vj >= Sij (3) 

The vectors u^ and v^ which are related by 

Avi = (Tsvd^i Ui (4) 

will be called normal signal and actuator modes, respectively, since they are eigenvectors of so called normal 
matrices, which are matrices which commute with their adjoints. The normal matrices are A^A in the case of 
the actuator vectors and AA^ in the case of the signal vectors. In the following the term normal mode without 
a specification 'actuator' or 'signal' will refer to an actuator mode. 

For a large number of rigcg of segments the normal modes with the smallest singular values become rather 
smooth. In the limit of an infinite number of segments over a finite area, that is in terms of the radii a of a 
segment (see figure 5) and R of the mirror, a/ R -^ 0, they should converge to functions which are continuously 
differentiable an infinite number of times. Let L denote the continuous operator which corresponds to the 
matrix A, and /i and gi the functions, called the analytical position and signal modes, to which the actuator 
modes Vi and the signal modes Ui, respectively, converge for a/R -^ 0. The position functions /i will also be 
referred to as analytical normal modes. The two now infinite sets {/i} and {gi} will have, similar to the sets 



{vi} and {ui}, the following two characteristics. Let the vector products be defined by 

1 /•+^ 1 '•+1 

'-D 



</.,/,>=^/ dxMx)f,{x) = - I dU^{Om) (5) 



for the one-dimensional problems with a length of the chain of segments of 2D and the normalised variable 
^ — x/D and by 

< f^Jj >= —r52 2T / d*^ / ^dr/,(r,(^)/j(r,(/?) = — j-/ d(^ / pdp f,{p,ip)fj{p,ip) (6) 

TTl^/t — r^j Jq Jr^ TT[L — pi) Jo Jpi 

for the two-dimensional problems, where r is the radial, (p the azimuthal coordinate, R the outer radius of 
a circular mirror, ri the radius of its inner hole, p — r/R the normalised radial variable, and pi — ri/R the 
normalised radius of the inner hole. The first characteristic is then that the two sets of functions are orthonormal, 

that is 

< ft, f] > = < 9^,9] >= ^^,j (7) 

The second characteristic is that the members fi and gi are related by 

L Ii == cFi gi , (8) 

where the (Xi, which from now on will be called analytical singular values, correspond to the singular values 
fsvd.i in the discrete SVD. 

The functions fi and gi can be derived from the two conditions (7) and (8). Introducing equation (8) into 
equation (7) gives 

^1 < 9z,9j > = < Lf„ Lfj > ^ a1 (5j,j (9) 

After integrating the vector product < Lfi, Lfj > by parts one gets [4] 

< L/„ Lf, > =< L"Lf,, /, > = < L^Lf,, f, > +^J/„ /,] , (10) 

where L^ is the adjoint of L, L^ the formal or Lagrange adjoint of L, and /3 [fi, fj] a boundary functional. If 
the boundary functional is zero equations (7), (9), and (10) lead to the linear differential equation 

L^Lf^d^f (11) 

Equations (10) and (11) are the fundamental equations for the derivation of the expressions for the analytical 
normal modes in the limit of infinitely fine segmentation. The boundary conditions which are required to derive 
specific infinite sets {/i} of analytical normal modes and the corresponding sets {di} of analytical singular values 
from the general solution of equation (11) are obtained from the requirement that the boundary functional 
vanishes. The form of the operator L^ and the boundary functional depend on the operator L. One of the 
problems is therefore to find the correct expression for the linear differential operator L which corresponds, in 
the limit of infinitely fine segmentation, to the matrix A. 



3 One-dimensional segmentation 

The method outlined in section 2 is first applied to the mathematically easier case of one-dimensional segmen- 
tation. In addition, some of the features appearing in this section are similar to the ones occuring in the more 
interesting examples of two-dimensional segmentation discussed in sections 4 and 5. 

3.1 SVD solutions 

Figure 1 shows two examples which are based on linear chains of segments. In figure la each segment can only 
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Figure 1: (a) : Segments with piston movements only, (b) : Segments linked at the interfaces. 

be moved in vertical direction. If the actuator positions are denoted by Zi^ the signals Si measured by vertical 
displacement sensors are given by 

Si = 2i+i - Zi (12) 

In figure lb the segments are linked and the actuators are at the links and at both ends of the chain. The 
signals Si are defined by 

Si^ Zi- 2zi+i + Zj+2 (13) 

For small angles Si is approximately proportional to the difference between the tilts of the segments. The 
matrices linking the actuator positions to the signals are given by 
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for the example (a) with rigig = riact — 1 and 
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for the example (b) with Usig = riact ~ 2. In both examples the number rznz of normal modes with non-zero 
singular values is equal to the number of signals. 

For large values of riact the lowest modes, which are the ones with the smallest singular values, approach 
smooth functions. Figure 2 shows for the example (b) with 100 segments the three modes with the lowest 
singular values. 




Figure 2: Example (b) ; Actuator values (circles) and sensor signals (crosses) for the three SVD modes with the lowest 
singular values. The solid lines are the mode shapes calculated analytically for an infinite number of segments. 



3.2 Analytical solutions for one-dimensional segmentation 

3.2.1 Example (a) 

In the limit of infinitely small segments the operator L which corresponds to the matrix Aa of example (a) is 
given by 

Then, with —D and D being the left and the right end of the chain of segments, equation (10) becomes 



1 / >2 X . r/ 1 X n+1 



^ ^ '■.niO]m)d^ 



^.m)]fAO 



For an arbitrary choice of a function fj from the set {fj} the boundary functional vanishes only if 



(18) 



The differential equation (11) in terms of the normalised variable ^ is then given by 

-^2$m = ^'f(o, (20) 

or, with the new normalised parameter a ~ Da which has the advantage that it is independent of the length 
ID of the chain and has got the unit 1, by 

d2 



'^.no = ^'f{o (21) 

The general solutions are either symmetric or antisymmetric : 

MO = c,cos«) (22) 

/a(e) - CaSin(aO (23) 

with Cs and Ca denoting arbitrary constants. Introducing these equations into the boundary conditions (19) gives 
the equations from which the infinite sets of normalised analytical singular values ai = Dai for the symmetric 
and the antisymmetric cases can be obtained. The ai are identical to the wavenumbers of the solution functions 
and given by 

<^^ - * I (24) 

Obviously, for the symmetric and antisymmetric cases the sets {/s/a,i} ^nd {i/s/a.i}, representing the analytical 
position and signal modes respectively, are orthogonal sets of functions. The following relationships, which follow 
from the equations (7), (9), and (17) will be required in section 3.3.1 : 

.2-2 _ n2 < Lfi{x),Lf^{x) > 



at = D' at = D 



< U{x)Ji{x) > 

+ 1 / J \ 2 .+1 



1 (^/«(^)) '^^/ J\ iMO)'d( (25) 



3.2.2 Example (b) 

In the limit of infinitely small segments the operator L which correspnds to the matrix A^, of example (b) is 
given by 









Taking the same steps as in section 3.2.1 one gets 
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The differential equation is therefore given by 
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/(e) = 



(26) 

(27) 
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(29) 



with a = D^a. The fourth-order differential operator in equation (29) can be factorised into two second-order 
differential operators : 

d^ d2 

(37o-<^)(377 + '^)/(0 = 



'Ae dC2 

With the definition 

the symmetric and antisymmetric solutions are, with as and a^ denoting additional coefficients, 

/s(e) = Cs[cos(A^) + asCosh(AC)] 
/a(e) = Ca[sin(AO+aaSinh(AO] 



(30) 

(31) 

(32) 
(33) 



Both sets of solutions require two boundary conditions which can be derived from the boundary functional (28). 
This must be zero for arbitrary functions taken from the sets of functions {fs.j} or {/a.j}. Since at least for 
some of them both the function and its derivative are non-zero at the edges the four conditions which guarantee 
that the boundary functional is zero for all combinations of i and j are 
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(34) 



(35) 



Introducing the general solutions (32) and (33) into equations (34) and (35) gives the equations from which 
the wavenumbers Ai and the coefficients as,i and a^^i can be calculated. The wavenumber Aj of a mode i is 
approximately given by 

A, == - (1 + 2i) (36) 



J(l + 2^) 



This approximation is very accurate for large i. For small i the differences are only of the order of a few percent. 
According to equation (31) the normalised analytical singular values Gi are the squares of the wavenumbers A^. 



3.3 Comparison between SVD-modes and analytical modes 
3.3.1 Example (a) 



One can derive an expression for the singular values crsvd,an,i expected from SVD based on the normalised 
analytical singular values (Ji. From the equations (3) and (4) one gets 

2 < Av,,Avj > < Avi, Avi > 



-^svd.i 



< Ui,Ui > 



< Vi,V, > 



< AV,:, AVi > , < Vi, V,; > 



(37) 
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^sig sig,i 

^act ^act.i 



where Csig.i and eact,i ai'e the r.in.s. values of the signals and the actuator positions, respectively. To derive an 
expression for esig,i let d be the length of a segment, rigeg the number of segments, and 2D = Uscgd the length of 
the full linear chain. If fi{x) is one of the analytical normal modes going through the centers of the segments 
and the Xj are the locations of the sensors, the mean square e^jg ^ of the signals is approximately given by 



g2 _ j2 
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For a very large number of segments the averaged sum can be replaced by the integral expressions in equation (17) 
with i = j, leading to the analytical mean square e'^- ^^ ^ of the signals : 

sig,an,< 2\Dj 7_i ^ dC 



4g,an,. = o 7^ ^ir^ de (39) 



Dividing equation (39) by the analytical mean square e^^^. ^^ ^ of the positions 
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4t,a„,.= o / /'(OdC (40) 



and using equation (25) gives 
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(41) 



Introducing equation (41) into equation (37) and taking the square root gives an expression not any more for 
the true singular values CTsvd,i but for the ones (Tsvd.an.i based on the analytical approximation with an infinite 
number of segments. The relationship between the normalised analytical singular value ai of the mode i and 
the corresponding one (Tsvd,an,i expected from SVD is then 



"■svd.an.i = \ ;=: Cj (42) 

V "act D 

For a chain with 100 segments figure 3a shows the singular values asvd.a,n,i according to equation (42) (crosses) 
and the true singular values crsvd,i calculated by SVD (circles). As one would expect, they only coincide for 
the lowest smooth modes. For the higher modes the analytical calculation overestimates the singular values. 
This can also be seen in figure 3b which shows the ratios of CTsvd.an.i to (Tsvd,i- Also for the lowest modes these 
ratios are not exactly equal to 1, but approximately equal to 1 + 0.5/nscg- For a large number of segments 
the shapes of the SVD solutions with the smallest singular values are in perfect agreement with the ones of 
the analytical solutions. Since the normal modes are used for the alignment of segmented mirrors, a common 
question is how the signal errors propagate to actuator errors, that is, which r.m.s. eact,n of the actuator errors 
is generated by noise in the signals with a given r.m.s. of esig,n- The r.m.s. eact,n can be calculated by a sum 
involving the Um non-zero singular values. Two r.m.s. values will be calculated. An analytical one denoted by 
fact, n. an, where the sum is formed with the singular values crsvd,an,i, and a true one denoted by eact.n,svd, where 
the sum is formed with the true singular values Usvd.i- In general, the mean square e^^.^. ^^ of the actuator errors 
based on analytical singular values is the sum of the mean squares e^^^. ^^ ^ of the individual modes. Using also 
equation (41) one gets 

e^ =Ve2 =Ve2. f —\ — (-431 

^act,an / ^ ^act,an,2 / ^ sig,an,z I J / _2 \ '^ / 
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Figure 3: Example (a), (a) : Singular values calculated analytically (crosses) and by SVD (circles), (b) : Ratios of 
these values for the lowest 50 modes. 

If the signal errors are random, each of the Usig signal modes contributes an equal amount egign.i = esig,!!/^/"' 
to the total r.m.s. esig,n of the signal noise. The analytically expected total mean square e^^.^ 



errors due to random signal noise is then, after replacing egig.an.i by ^s 
expression in equation (24), given by 



Csig,] 

of the actuator 
Csig.n/^/'T-sig and (7i by the explicit 



e^ 1- 

sig,n 2 



2D 



— e4 

SIS - 1 *" 
^ I— 1 



(44) 



For a large number of segments rigig can be replaced by 2D/d. Furthermore, because of its fast convergence, 
the sum from 1 to Unz is, for large riaz, approximately equal to the sum from 1 to infinity, with a value of 7r^/6. 
All this finally gives 

T f2D 



*^act,n,an ^^ ^sig,n \ r. \ i V ^/ 

V 6 V a 

For a given r.m.s. of the signal noise the r.m.s. of the actuator error therefore scales, as one would expect, with 

the square root of the number of segments. Equation (44) can also be written as 
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Similarly, the true r.m.s. eact.n,svd calculated with the true SVD singular values asvd,i is given by 

2 _!_\^ L 

sig,n „_^ Z^ ^2 ^ 



^act.n.svd 



^act .^^ ^svd,^ 



(46) 



(47) 



The result from equation (47) is effectively identical to the one from equation (45). The same result is also 
obtained by simple statistical considerations, that is by summing up random piston movements along a line of 
segments and calculating the r.m.s. after subtracting the average. 



3.3.2 Example (b) 



The expression for asvd.an.i can be derived in the same way as in section 3.3.1. One gets for a specific mode i 
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with (Ti = D^&i. Here egig.an.i is the r.m.s. of the angular signal error in radians. The singular value tTgvd,: 
expected from SVD is then related to the corresponding analytical singular value ai by 



^svd,i 



JT-act 



D2 



(50) 



For a chain with 100 segments figure 4a shows the square roots of the singular values asvd.anA (crosses), that 
is the wavenumbers corresponding to (Tsvd,an,i according to equation (31), and of the true ones crsvd,i (circles). 
Figure 4a for the wavenumbers is very similar to the corresponding figure 3a for the singular values of example 
(a). Also the figure 4b for the ratios of the square roots of (7svd,an,i to the square roots of asvd.i is effectively 
identical to the corresponding figure 3b of example (a) for the ratios of the singular values. The total mean 
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Figure 4: Example (b). (a) : Wavenumbers calculated analytically (crosses) and by SVD (circles), (b) : Ratios of these 
values for the lowest 50 modes. 



square of the actuator errors is given by, using also equation (36) for the analytical singular values. 
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(51) 



Next, risig is replaced by 2D/d. Furthermore, because of its fast convergence, the sum from 1 to Unz is, for large 
rinz, approximately equal to the sum from 1 to infinity, with a value of 7r*/120 — 1. All this finally gives 



eact.n.an « 0.0982 D \j^ £, 

V « 



(52) 



For a chain of 100 segments the difference between the true value eact,n,svd obtained from equation (47) applied 
to example (b) and eact,n.an from equation (52) is less than 1%. The solid lines in figure 2 show the analytical 
solutions with the three smallest singular values. Clearly, they are in perfect agreement with the SVD solutions. 



4 Analytical solutions for two-dimensional segmentation 



4.1 General remarks 



There are several combinations of signals and actuator positions which could be studied. In this chapter we 
discuss mainly two examples which are of particular interest for the alignment and phasing of segmented mirrors. 



We assume that the segments are hexagonal and the signals are measured at locations along the intersegment 
edges as shown in figure 5. In both examples the signals are the relative vertical displacements between adjacent 
segments at the locations of the sensors. Such signals can be obtained by optical measurements or by position 
sensors. In principle, the signals can also include the relative tilts between adjacent segments, defined here as 




Figure 5: Locations of the sensors at the intersegment edges of hexagonal segments. The filled circles in the segments 
mark the positions of the actuators. 

a change in the slope of the segments perpendicular to the edges at the locations of the sensors. Usually, the 
signals generated by relative tilts are much smaller than signals generated by relative vertical displacements [6] . 
Except for one particular mode, the dcfocus mode, which is discussed in section 5.4, relative tilts at the sensors 
are therefore of less significance. 

The two examples differ in the number of the degrees of freedom for the movements of the segments. 

• Piston movements 

In this case, which is the two-dimensional equivalent of the one-dimensional example (a) presented in 
section 3.2.1, the segments can only move in piston and the positions of the three actuators of one 
segment arc therefore all identical. The signals along one intersegment edge are all the same and it would 
be sufficient to consider just one signal along one edge. This example describes a situation where the 
segments have already been perfectly aligned in tilt and the remaining phasing errors are corrected by 
pure piston movements of the segments, a procedure which is used for the phasing of the segmented 
primary mirror of the Keck telescope based on optical measurements. 

• Piston, tip and tilt movements 

In this case, which is to some extent related to the one-dimensional example (b) presented in section 3.2.2, 
there are three degrees of freedom for the movements of the segments, namely piston, tip and tilt. It 
describes situations where a full correction of the misalignments including phasing is based on relative 
vertical displacement signals at the locations shown in figure 5. It is used for the fast stabilization of the 
alignment in the Keck telescope with signals obtained by position sensors. 



The sensors at the edges of the hexagonal segments detect signals in the three directions 9i = —60°, 6*2 = 0° 
and ^3 = 60° with respect to the horizontal axis in figure 5. The derivative of a function /(r, ip) in one of the 
directions 6k, k — 1,2, 3, expressed as derivatives with respect to the radial and azimuthal coordinates r and (p 
of a cylindrical coordinate system, is given by 
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(53) 
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In the limit of infinitesimally small segments the lowest SVD modes always follow rotational symmetries n. 
They can therefore be written as a product of a radial and an azimuthal component : 



f{r,y^) = fn{r) cos{nip) 



(54) 



The differential operators for curvatures are second derivatives based on the operator in equation (53). Here 
we use the following approximation. Strictly speaking, the curvature K along a direction s is defined by 
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i+'d7 



-1 3/2 



(55) 



Since the derivative in the denominator is small compared to unity, we will replace the denominator by 1. This 
is equivalent to approximating the spherical surface by a paraboloid with a radial dependence on p^. 



4.2 Piston movements 

The two-dimensional case which is equivalent to the one-dimensional example in section 3.2.1 is a configuration 
where the segments can only perform piston movements and the sensors measure only the relative vertical 
displacements. In this case the differential operator is given by the operator Tk in equation (53). The starting 
point for the derivation of the analytical normal modes is a vector product averaged over the three directions 
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After introducing equation (54), integrating over the azimuth angle (p, and summing over k equation (56) 
becomes 
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Following the procedure outlined in section 2 one gets 
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where Lb,™ is the operator appearing in the Bcssel differential equation : 

_ d2 Id n2 

LB,n — -J— J H "5 J 
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The general solution of 

with a related to a in equation (11) by a — \/2Ra is a sum of the Bessel functions Jn{ap) and Yn{(yp) with 
coefficients cj and cy : 

/n(p) = cjJn{ap) + cyYniap) (62) 

For arbitrary functions fn,j taken from the set of general solutions (62) the boundary functional (59) vanishes 
if the expression in the square brackets is zero for both p = pi and p = 1, that is if the derivatives vanish 
at both the inner and outer edges of the mirror. Introducing the general solution (62) into the two boundary 
conditions gives a set of two linear homogeneous equation in cj and cy from which the infinite set of normalised 
analytical singular values an,i and the corresponding ratios {cy /c])n,i can be obtained. As in the corresponding 
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one-dimensional example discussed in section 3.2.1 the singular values an,i are identical to the wavenumbers of 
the functions, in this case two Bessel functions, contributing to the solution (62). 

For the rotational symmetries 0, 1, and 2 figure 6 shows, for pi = 0.1, the four analytical normal modes 
with the smallest singular values an,i- Like the Zernike polynomials the modes can be classified according to 
their rotational symmetries and, within each rotational symmetry, according to their number of nodes, that is 
zeros along the radial coordinate. The only mode, which is, compared with the set of Zernike modes, missing in 
the set {fn,i\ is the lowest mode of rotational symmetry with no nodes. This reflects the fact that a uniform 
piston movement of all segments does not generate relative vertical displacements at intersegment edges. 

The major difference to the Zernike polynomials is that, due to the boundary conditions following from 
equation (59), the first derivatives are, as in the one-dimensional example in section 3.2.1, zero at the edges. In 




Figure 6: Lowest analytical normal modes of an annular mirror with pi = 0.1 for pure piston movements for the 
rotational symmetries 0, 1, and 2. 

the case of a completely filled mirror, that is pi — 0, the function Yn{<Tp), which diverges at p = 0, drops out of 
the general solution, which is then given by 



fn{p) ^C,TJn{(7p) 



(63) 



In this case the expression inside the square brackets of the boundary functional (59) is automatically zero 
for p = 0. This leaves therefore only one boundary condition which is sufficient to determine the infinite set 
of normalised analytical singular values (T„.i which are also identical to the wavenumbers. The functions look 
similar to the ones in figure 6, but now, for all symmetries n > 0, they are zero at the centre of the mirror. 



4.3 Piston, tip and tilt movements 

4.3.1 Differential operator 

As in section 4.2 the sensors detect only relative vertical displacements, but the segments can now be controlled 
in piston, tip and tilt. Let a smooth function over hexagonal areas be approximated by planes. For a spherical 
surface this can be done without generating relative axial displacements. The sensors therefore measure locally 
only the difference to a surface with, also locally, a constant curvature in all directions. This is the same as the 
torsion of the surface which is given by, using also equation (54), 



Lkf{r,ip) 



—^ [sn\{2{ip - Ok)) cos{n(p) Ln,ifn{p) + "2 cos{2{ip ~ 9k)) sin{nip) Ln,2fn{p)] (64) 
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with the two operators Ln,i and Ln.2 defined by 



L 



(P 1 d 
dp^ p dp 
n d n 



n.l — — 



pdp 



(65) 
(66) 



4.3.2 Derivation of the analytical normal modes 

The procedure outlined in section 2 starts off with the following vector product averaged over the three direction 
Ofc : 

3 -, -, 3 , 

-kMp,f))iLkfjip,fj) (67) 



1 1 1 /■ /"^ 

- ^ < Lkfi{r,(p),Lkfj{r,(p) > = -— jT ^ h^^ P<ip{Lk 

•i fe=i •^ ^y^ Pi) k^iJ J PI 



After introducing equation (54), integrating over the azimuth angle (^3, and summing over the three directions 
denoted by k one gets 



1 ^ 

- ^ < Lkji{r, (f), Lkfj{r, f)> = 



fc=i 



1 1 + Sn,a 



Pi -J PI 

4(i„,2/«,j(p))(i«,2/«j(p))] 



pdp [{Ln,lfns{p)){Ln,lfn,j{p)) + 



A straightforward calculation shows that 



The boundary functional is given by 

1 1 + 5n.O 



^^^"8i?^^^' 



(68) 
(69) 



PT[Jn,iJn 



Pi 



p-^fnAP) 



dp 



UAp) -^ — fnAp) 



1 d 



dp 



InAP) 



Pi 



-p-T^fnAp) - T^fnAp) + (1 + 3n )-—fnAp) ~ 4 — /„,,(p) fn,]{p) 



P^ . 



dp-^ dp^ p dp p 

With a related to a in equation (11) by a = \/8R'^a the fourth-order differential operator appearing in 

Li,n.fn{p) =CrVn(p) 

can be written as a product of two second-order differential operators. Equation (71) then becomes 

{Lb,„ + cr) [Lb^u - 'y).fn{p) = 



(70) 



(71) 



(72) 



With the same definition (31) of A as in section 3.2.2 the general solution of (72) is a sum of the Bessel functions 
Jn{Xp), Yn(Xp), In{Xp), and Kn{Xp) with coefficients cj, cy, cy and ck : 



f{p) ^ Cj J„(A/9) + CyYniXp) + Ci/„(Ap) + CKKn{Xp) 



(73) 



For an arbitrary choice of the functions /„j the boundary functional (70) is zero if the expressions in both square 
brackets vanish for p = pi and p = 1. Introducing (73) into the four boundary conditions gives a system of 
four linear homogeneous equations in cj, cy, cj, and ck from which the infinite set of normalised wavenumbers 
Xn,i and the corresponding ratios of cy, ci, and ck to cj can be obtained. Similarly to the corresponding 
one-dimensional example the singular values an,i are the squares of the wavenumbers Xn,i- 

For the rotational symmetries 0, 1, and 2 and pi =0.1 figure 7 shows the four analytical normal modes with 
the lowest singular values an,i- Using the classification of the modes introduced in section 4.2 three modes are 
missing compared with the set of Zernike modes. Two of them are the two lowest modes of rotational symmetry 
without and with one node, corresponding to Zernike piston and defocus, and the third one is the lowest 
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mode of rotational symmetry 1 corresponding to Zcrnikc tilt. All three modes can be generated without relative 
vertical displacements at intersegment edges. 

The functions are somewhat similar to the corresponding Zernike polynomials, but should better be compared 
with the elastic minimum energy modes of circular plates (see section 7). For the case of a full mirror, that is 




Figure 7: Lowest analytical normal modes of an annular mirror with pi = 0.1 for piston, tip and tilt movements for the 
rotational symmetries 0, 1, and 2. 

pi = 0, the Bessel functions y„(Ap) and Kn{Xp), which diverge at p = 0, drop out of the general solution, which 
is then given by 

/-c,TJ„(Ap)+ci/„(Ap) (74) 

The expressions in the two square brackets of the boundary functional (70) are automatically zero for p ~ 0. This 
leaves two boundary conditions which are sufficient to calculate the wavenumbers Xn,i and the corresponding 
ratios {ci/c])n,i- The functions look similar to the ones in figure 7, but now, for all symmetries n > 0, they are 
zero at the centre of the mirror. 



5 Comparison with SVD results 

5.1 SVD programs 

All SVD calculations were done with computer programs supplied by G. Chanan [5]. They were developed for 
and applied to the control of the segmented primary mirror of the Keck telescope. Most of the data used in this 
paper were obtained for the segmented mirror shown in figure 8. In the following, the radius of the segments 
is, as shown in figure 5, denoted by a and the distance between the sensors along one intersegment edge by b. 
For the numerical SVD-calculations we used the parameters a — 0.5 m, b = 0.3078 m, and R = 15.05 m. The 
segments in figure 8 are the ones with the centres inside the range of normalised radii p of 0.30 < p < 1.0. This 
mirror has got 1002 segments and 5688 sensors. The radius of the circle Tact defined by the actuator locations 
shown in figure 5 was ract — 0.52275 a. Calculations were also done for a geometry similar to the one shown in 
figure 8, but without a central hole. This mirror has got 1099 segments and 6348 sensors. 



5.2 Piston movements 



5.2.1 Mode shapes 

The number of SVD modes with non-zero singular values is rinz — rig. 

equal to zero, which corresponds to the missing analytical mode in section 4.2, is related to a uniform piston 



1. One mode with a singular value 
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Figure 8: Segmentation used for the SVD calculations. 

movement of all segm.ents. For a few low order modes figure 9 shows the traces along a diameter for a full 
and an annular mirror. Especially for the full mirror the agreement between the shapes of the normal modes 
obtained with SVD (circles) and analytically (solid lines) in section 4.2 is excellent. The deviations close to 
the inner hole of the annular mirror are to be expected since this edge is less accurately defined because of the 
smaller number of segments. 

5.2.2 Singular values 



For pure piston movements a sensor signal is proportional to the difference between the values of an analytical 
normal mode fn,i at adjacent segment centres. To the first approximation, this is proportional to the slopes of 
fn.i midway between two centres along the line connecting the two centres times the distance y'Sa between the 
two centres. The mean square e^; „ ^ of the sensor signals of the mode (n, i) is therefore 



2 1 * 



rj^Z 



l-Slg 






(75) 



where Tnjj is the slope of the mode i of rotational symmetry n perpendicular to the edge at the location of 
the sensor j. For a very large number of segments the averaged sum in equation (75) can be replaced by the 
integral expressions in equation (57) with i = j : 



3a2 1 



2i?2 1 _ 



pI 



pdp 



dp 



In,i{p) 



After dividing by the mean square e^^ a.nn i ^f the actuator positions 



rfn.tip) 



(76) 



act,n,i 



pI 



P^PfnAP) 



Pi 



(77) 
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Figure 9: Normal modes calculated by SVD (circles) and analytical normal modes (solid lines) for pure piston movements. 
The first index along a curve denotes the rotational symmetry and the second one the order within the symmetry, (a) : 
Full mirror, (b) : Annular mirror with pi — 0.3. 



one gets 



sig,an,n,z 



act,an,n,i 



3a^ 
2i?2 



2 _ o 2 ~2 



(78) 



where an,i = ylRdn^i are the normahsed analytical singular values. Similarly to section 3.3.1 one finally gets 
for the relationship between a singular value asvd,an,n,2 expected analytically from SVD and the corresponding 



f^svd,; 



O /tgicr CI 

2nact R 



(79) 



Since the ratio ngig/nact is effectively constant for highly segmented mirrors, the singular values for a specific 
mode (n, i) expected from SVD scale with a/R. Figure 10a shows for the lowest 25 normal modes the true 
singular values asvd,n,i given by SVD (circles) on the one hand and the ones CTsvd, an, n,i (crosses) according 
to equation (79) on the other hand. The agreement is quite satisfactory. Figure 10b, displaying the ratios of 
Csvd, an, n,i to crsvd,ri,i, shows that for the higher modes the analytical calculation overstimates the singular values. 

5.2.3 Error propagation 

Starting from equation (78) one can derive how the total r.m.s. eact,n,an of the actuator errors generated by a 
given r.m.s. esig,n of the signal errors scales with the number of segments. The procedure is the same as the 
one used in section 3.3.1. One gets 



2i?2 1 



1 



£act,n,an ^sig,n o„2 „ . 2^ ^2 ' 



(80) 



where the sum runs over the n 

the total number of segments is approximately given by 



nz "'scg 



1 non-zero analytical singular values. For a large number of segments 



2tt (R 
''"' "* V3 U 



(81) 



Neglecting the missing sensors at the edge of the mirror, the number of sensors is approximately given by 



Tisig ''^ Tlscg 



(82) 
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Figure 10: Singular values calculated by SVD (circles) and analytically (crosses) for pure piston movements, (a) : 
Lowest 25 modes, (b) : Ratios of analytically expected singular values to the ones calculated by SVD for the lowest 230 
modes. 



Introducing all this into equation (80) one gets 



e:.„^^-L (83) 



'^'^^" 2^V3 



For the higher modes the analytical singular values increase approximately linearly both with the symmetry n 
and the order i, in the latter case, as one should expect, in steps of tt. With such linear dependencies on n and 
i the sum in equation (83) does not converge if the number of segments goes to infinity. But it increases only 
very slowly, namely approximately linearly with the logarithm of nnz. Therefore, eact,n.an scales approximately 
with log{R/a). 

For the full mirror desribed in section 5.1 the sum in equation (83) is approximately equal to 2.4. This gives 
for the r.m.s. eact,n,an of the actuator errors 

eact,„,an ~ 0.470 e^ig^n (84) 

This is in good agreement with the true r.m.s. eact,n,svd — 0.504 egig.n of the actuator errors calculated with 
the true singular values <Jsvd,n,i- Similar values for eact,n,svd are given in table 4 of reference [6]. The true 
value Cact.n.svd IS slightly larger than eact,n,an since, first, the singular values <Jsvd,a.n,n,i are larger than the true 
ones (Tsvd. n,i and, second, the number of signals is smaller than the value of 6 rigeg used in the derivation of 
equation (84). 

5.3 Piston, tip and tilt movements 

5.3.1 Mode shapes 

The number of non-zero singular values is Unz — J^act — 4. Four normal modes calculated by SVD with singular 
values zero are the piston mode, the defocus mode and the two tilt modes. They correspond to the three missing 
analytical normal modes in section 4.3.2. Figure 11 shows for a few of the lowest modes the SVD solutions 
(circles) and the corresponding analytical solutions (solid lines) for a full and an annular mirror. Especially for 
the full mirror the agreement is excellent. 
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Figure 11: Normal modes calculated by SVD (circles) and analytical normal modes (solid lines) for piston, tip and tilt 
movements for the lowest two modes of the rotational symmetries 1 and 2. (a) : Full mirror, (b) : Ammlar mirror with 
pi = 0.3. 

5.3.2 Singular values 



The sensor signals are proportional to the difference in the slope of an analytical normal mode /„.i between 
adjacent segment centres. If only the torsion contributes to a sensor signal one has to take the slope perpendicular 
to the line connecting the two centres. The difference in the slope is then, to the first approximation, proportional 
to the torsion of ]n,i at the location midway between the two centres, the distance v3a between the two centres, 
and the distance from the midway point to the sensor, that is 6/2. The mean square Cgj „ ^ of the sensor signals 
for the mode (n, i) is therefore 



\fZah 



1 



(85) 



where Tn,i,j is the torsion of the mode (n, i) at the location of the sensor j. For a very large number of segments 
the averaged sum in equation (85) can be replaced by the integral expressions in equation (68) with i = j : 



^sig,an,n,z 



32i?4 T 



pI 



pdp[{L„^lfn4p)f + 4:(Ln,2.fn4p)? 



(86) 



After dividing, as in section 5.2.2, by the mean square of the actuator positions, the ratio of the mean square 
of the sensors signals to the mean square of the actuator positions for a given mode (n, i) becomes 



sig,an,n,z 



act, an, n, 2 



80252 2 



80252 ^ 



(87) 



where cr„.i — ^/%E?dn,i are the normalised analytical singular values. Similarly to section 3.2.1 one finally gets 



for the relationship between the singular value a^vd,: 
the corresponding normalised analytical one cr„^i 



of the mode (n, i) expected analytically from SVD and 



C^svd,; 



3 nsitr ah 



32 nact i?2 



Since the ratio nsig/riact is effectively constant for highly segmented mirrors, the singular values expected from 
SVD scale with ab/R^. Figure 12a shows for the lowest 25 normal modes the true singular values <Tsvd,n,i given 



by SVD (circles) on the one hand and the ones CTsvd,^ 



(crosses) according to equation (88) on the other 



hand. The agreement is quite satisfactory. Figure 12b, displaying the ratios of (Jsvd,i 
for the higher modes the analytical calculation understimates the singular values. 



to (Tgvd n i I shows that 
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Figure 12: Singular values calculated with SVD (circles) and analytically (crosses) for piston, tip and tilt movements. 
(a) : Lowest 25 modes, (b) : Ratios of expected singular values to the ones calculated by SVD for the lowest 230 modes. 

5.3.3 Error propagation 

Starting from equation (87) one can derive an expression for the total r.m.s. fact, n, an of the actuator errors 
based on the analytical singular values. The procedure is again the same as the one used in section 3.3.1. One 
gets 



3a262 



where the sum runs over the rinz = riact — 4 non-zero analytical singular values. For the higher modes the square 
roots of the analytical singular values increase approximately linearly both with the symmetry n and the order 
i, in the latter case, as one should expect, in steps of tt. The sum in equation (89) converges rapidly for a 
large number of segments. For the full mirror described in section 5.1 the limit value is approximately 0.065. 
Introducing this and the approximate expressions (81) and (82) into equation (89) one finally gets for the r.m.s. 
of the actuator errors 

D 

^act,n,an ^^ U.oi — (^sig,n W^j 

The true r.m.s. cract,n,svd of the actuator errors is effectively the same as the one calculated from equation (90). 
For the currently envisaged extremely large telescopes the ratio R/b is approximately equal to 50. The r.m.s. 
of the actuator errors is then approximately fifteen times larger than the r.m.s. of the signal errors. This is in 
good agreement with the results of SVD calculations for a thirty meter telescope given in table 1 of reference 

[6]. 

5.4 Piston, tip and tilt generated by relative tilt signals 

Some sensors detect in addition to the relative vertical displacements also the relative tilts between adjacent 
segments. One could use the same formalism as in the other examples to derive the differential equations and 
the boundary conditions. Unfortunately, the fourth-order differential operator cannot easily be factorised into 
two second-order differential operators as it could be done in equations (71) and (72). But usually the signals 
due to the relative tilts perpendicular to the edges are much smaller than the signals due to relative vertical 
displacements [6]. Most modes will then hardly be affected by the tilt signals. There is only one additional 
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mode of rotational symmetry zero with a non-zero singular value which converges, in the limit of the signals 
due to relative tilts going to zero, to the pure defocus mode with a singular value equal to zero. For small, but 
still non-zero tilt signals the singular value of this mode can be estimated. It is related to the integral of the 
mean curvature over the mirror for a mode shape proportional to pure defocus. The operator which measures 
the curvature along a direction 9^. is given by 



^'■^'-•-'•^ ^ i^^"™'"^> 



2tQ .9 fnip) , .. 2//J n1 (dfn{p) 



cos\e,)^^jfp. + sin^(0fe)- :^^^ - - Up) 

Op'^ p \ op p 



+ sin(n^) sin(20fe) n (^^ - - Up)) | (91) 

One then obtains, after introducing equation (54), integrating over the azimuth angle ip, and averaging over the 
three directions denoted by k, for the special case of rotational symmetry and a full mirror 

'd\fo4p)V , 3 fdfoAp)V , 2d^fo4p)dfo4p)' 



2 3 1 /"^ 

- ^ < Lfe.A(r, (/3),Lfc/j(r, (^) >= ^^ / pdp 



3 



dp^ J P^ \ dp J p dp^ dp 

(92) 
If one assumes that the lowest mode /o,2 (/o,i would denote the piston mode) is equal to pure Zcrnike defocus 
cV3(2p^ — 1), where c is the r.m.s. of this mode, one gets 

1 J^ ^ 48 

- ^ < LkfoAr, (/?), Lfc.fo,2(r, (/?) >= c^ — (93) 

k=l 

The sensor signal is proportional to the difference in tilt between adjacent segment centres. If only the curvature 
contributes to the sensor signal one has to take the tilt component parallel to the line connecting the two centres. 
The difference in tilt is then, to the first approximation, proportional to the curvature at the location midway 
between the two centres, the distance \/3a between the centres and the lever g which converts the angular 
difference into a differential displacement [6] . The mean square e'^- q 2 "^^ ^^'^ sensor signals is therefore 

"sig 

4,0,2 = i^agf -— E ^0.2,, , (94) 

'^'g j=l 

where To^2,j is the curvature of the defocus mode at the location of the sensor j. In the limit of small segments 
the averaged sum in equation (94) can be replaced by the integral expression in equation (93) : 

4g,a„,o,2 - 144 c^ (^)' (95) 

After dividing, as in section 5.2.2, by the mean square c^ of the position function, the ratio of the r.m.s. of the 
sensor signals to the r.m.s. of the actuator positions for the defocus mode becomes 

£sig,an,0,2 ^ ^2 ^ (96) 

The singular value expected from SVD is then finally given by 



CTsvd.an.0,2 — -L^ l/ 777 [^ ' ) 

V "^act R^ 

For the full mirror described in section 5.1 with g ~ 23.53 mm one gets from equation (97) (Tsvd,an,o,2 = 0.0007206. 
This is in good agreement with the result of a^vd 0.2 ~ 0.0007064 obtained from a SVD calculation. 



6 Relationship to modes minimizing the signals measured by the 
sensors 

For arbitrary position functions the mean square J^ of the signals is given /(^) by 

J=^ JiLfiOfd^ (98) 
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for the two one-dimensional examples with the operators L defined in equations (16) and (26), and for arbitrary 
position functions fn{p) by 



Pi J PI \ dp P 

^ = ^V^ / /'dp((iM/n(p))'+4(ip,2/n(p))') (100) 

for the two-dimensional examples in sections 4.2 and 4.3. On the assumption that the mean square C of the 
position functions, defined by 

C = \ r fiOd^ (101) 



^ - TH^^^^ J^/<ip[(^/nip)r + ^.fM] (99) 

1 1 + '^n.O 
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n.O / J #2 



C = -r-^ / Pdp/^(p) (102) 

for the one- and two-dimensional cases, is constant, the modes which minimise the r.m.s. of the signals can be 
calculated by variational methods. Denoting the variation by S one has to solve 

6{J - CC) - (103) 

where C is a free parameter which can be interpreted the energy per unit of the r.m.s. of C This gives, with 
m — and rj — ^ for the one-dimensional and to = 1, 77 = p, and f = fn for the two-dimensional examples 

^ - 2Cry™/ = (104) 

or explicitly, with a prime denoting a partial derivative with respect to 77, 

For all examples given in this paper the differential equations derived from equation (105) are identical to the 
ones derived with the method outlined in section 2. The analytical normal modes can therefore be regarded as a 
special class of the modes which minimise the r.m.s. of the signals for a specified r.m.s. of the position function 
/. Their special characteristic is that the position modes and corresponding signal modes form orthogonal sets, 
which defines the boundary conditions. 



7 Relationship of analytical normal modes to elastic modes 

7.1 Reasons for a comparison 

In future extremely large telescopes segmented mirrors will be used together with large fiexible meniscus mirrors. 
Whereas the aberrations of the segmented mirrors are best described by the analytical normal modes derived in 
this paper, the deformations of the thin meniscus mirros are best described by minimum elastic energy modes 
[3]. For the two examples in sects. 3.2.2 and 4.3 both types of modes are closely related and, as will be shown 
in this section, often very similar. This offers the possibility to correct aberrations, which are generated by 
segmented mirrors, by thin meniscus mirrors and vice versa. 

7.2 One-dimesional segmentation 

The differential equation (30) for example (b) in section 3.2.2 is identical to the one for minimum elastic energy 
modes of a bar. Since the boundary conditions (34) and (35) are also the same as the boundary conditions for 
a bar which is free at both ends, the analytical normal modes of example (b) are identical to the elastic modes 
of a free bar. 
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Table 1: Residual r.m.s. after fitting a set of elastic modes of rotational symmetry 2 to the four lowest analytical normal 
modes of rotational symmetry 2. 
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0.0056 
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Normal 
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0.9969 


0.0053 


0.0007 
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0.0006 
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0.9998 


0.9998 


0.0124 


0.0066 


0.0042 
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1.0000 


1.0000 


0.9999 


0.0156 


0.0081 



7.3 Two-dimensional segmentation 

It is shown in appendix A that the orthogonal set of minimum elastic energy modes of a thin circular plate can 
be derived with a slightly more general method than the one introduced in section 2. For the set of elastic modes 
the corresponding set of curvature tensors has to be orthogonal to the corresponding set of moment tensors. 
This leads to the same differential equation (72) as for the case of piston, tip and tilt movements in section 4.3.2. 
Any differences in the shapes must then be related to the boundary conditions. For higher order modes with 
large wavenumbers the terms with highest derivatives in the boundary conditions dominate the other terms. 
If only the terms with the highest derivatives are retained, the boundary conditions for the analytical normal 
modes in equation (70) and the elastic modes in equations (117) and (118) below of a free plate are identical : 



d3 
Jp3MP) 



p=i 



p=i 



(106) 
(107) 



Therefore, within each rotational symmetry, the higher analytical normal modes should be very similar to the 
corresponding elastic modes of a circular plate. This similarity is obvious from figure 13 which compares the 
shapes of the lowest analytical normal modes of the rotational symmetries and 2 with the shapes of the 
corresponding elastic modes. 




Figure 13: Comparison of the lowest analytical normal modes (solid lines) and corresponding clastic modes (dashed 
lines) for the rotational symmetries and 2. 



Table 1 shows, for an initial r.m.s. of an analytical normal mode of 1, the residual r.m.s. after fitting a given 
number of the lowest elastic modes of a full circular plate to the lowest four analytical normal modes of rotational 
symmetry 2 of a full mirror. The difference between the lowest analytical normal mode and the corresponding 
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lowest elastic mode is of the order of 8% and one needs to fit the lowest three elastic modes to push the residual 
r.m.s. to a level of 1%. The similarity of the higher order modes is expressed by the fact that the residual r.m.s. 
is effectively equal to 1 if the lowest z — 1 elastic modes are fitted to i-th analytical normal mode, but drops 
sharply to approximately 1% after the i-th elastic mode is included in the fit. 

It is interesting to compare the probabilities of the occurence of normal modes in a segmented mirror on 
the one hand with the probabilities of the occurence of elastic modes in a monolithic mirror on the other hand. 
One has to choose the normal modes obtained from the more general case of piston, tip and tilt movements 
presented in section 4.3 because of their similarity to the elastic minimum energy modes. 

The source for wavefront errors in the form of normal modes are the signals at intersegment edges. For 
white noise the r.m.s. values esig,i of all signal modes are identical. Since the singular value is proportional to 
1/A| the r.m.s. enor,i of a normal mode is related to Csig by 

e„or,^ CX ^ (108) 

\ 

The source for wavefront errors in the form of elastic modes are errors in the support forces of the large 
monolithic mirrors. For a large number of supports these force errors can be replaced by pressure fields. Since 
the elastic modes form an orthogonal set of functions and are proportional to the pressure fields generating 
them, the elastic modes and pressure functions can be regarded as a pair of normal modes similar to the pair 
of analytical signal and position modes in segmented mirrors. But, instead of being proportional to 1/Af the 
r.m.s. eda.i of an elastic mode is, for a given r.m.s. Cpr of the pressure field, proportional to 1/A^ : 



A^ 



eeUa cx ^ (109) 



A comparison of the equations (108) and (109) shows that the coefficients of the elastic modes generated by 
random force errors decline much faster with the order of the mode than the coefficients of the normal modes of 
a segmented mirror generated by signal errors. The difference in the behaviour is, at a first glance, surprising, 
since the normal modes of a segmented mirror are similar to the elastic modes, that is the relationship between 
the position modes and the fields of the second derivatives of the position modes is the same. But the difference 
is that for the normal mirror modes the source is already a kind of curvature field, whereas for the elastic modes 
there is another additional relationship involving a factor Xf between the source, which is the pressure field, 
and the curvature field. 

Appendix A : Derivation of elastic modes of a circular plate 

For the derivation of the elastic modes of a free circular plate one has to generalise the method introduced in 
section 2. Instead of equation (9) 

< L/„ L/, > ex (5,,j (110) 

we require now that 

< L/„ TV/, > ex ,5,,, (111) 

with two possibly different differential operators L and A^. Partial integration of the left hand side leads to an 
equation which is equivalent to equation (10) : 

< L/„ Nf, > - < N^L,U, f, > +/3l,n[/- /j] (112) 

If the boundary functional /3 ^[fi,fj] is zero, equations (7), (111), and (112) lead to the linear differential 
equation 

N^'Lf, = erf /, (113) 

For elastic modes the operators are related to the strain and moment tensors. With a prime and an asterix 
denoting a partial derivative with respect to the radial variable r and the azimuthal variable (j), respectively, 
the components of the strain tensor are defined by 

Xr--/", X0--7-7F' Xr0 = -^ + ^ (114) 
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and the components of the moment tensor by 

M, = D{xr + i^X<t>), M4, = D{x4. + vxr), Mr^ = D{1 - v)xr4, , (115) 

where D = Kh^ /12 and K = _E/i/(l — z^^). E is the modulus of elasticity, h the mirror thickness, and ly Poisson's 
ratio. Equation (112) then becomes 

<Lfi,Nfj> =^ dip dr {Mr^iXrj + M^^iX4,,j + 2Mr^^iXr4,,j) 

Jo Jri 

fr2 

dr [{Xr,iXrj + X4>M4>-J + 2Xr0,*Xr0,i) 
'ri 

+ l^ {X4>,iXr,j + Xr,iX<t>,3 ~ '^Xrit>,iXr<t>,j)] 

oc 5ij (116) 

Introducing equation (54) and partial integration lead to the differential equation (72) and a boundary func- 
tional. For an arbitrary choice of the functions fn,j the boundary functional is always zero if the following 
boundary conditions are fulfilled : 

f^{r^) + J-l^-n^J-^ = (117) 

r:\r,)+^-!^ + {l + n\2-v))^-^+nHi-v)^-^ = (118) 

for both Vx — ri and r^ = ^2- These four boundary conditions are identical to those for a free thin annular 
circular plate [3] . 
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